Multiplicity and concentration of solutions for a fractional $ p $-Kirchhoff type equation

This paper is concerned with the following fractional $ p $-Kirchhoff equation \begin{document}$ \begin{eqnarray*} \left\{ \begin{array}{ll} \varepsilon ^{sp}M\left( {\varepsilon ^{sp - N}}\iint_{\mathbb{R}^{2N}}\frac{{{{\left| {u(x) u(y)} \right|}^p}}}{{{{\left| {x y} \right|}^{N + sp}}}}dxdy\right)(-\Delta)_p^su V(x){u^{p 1}} = {u^{p_s^* 1}}+f(u), \\ u>0, \ \mbox{in}\ {\mathbb{R}^N}, \end{array} \right. \end{eqnarray*} $\end{document} where \varepsilon>0 a parameter, M(t) a+bt^{\theta-1} a>0 $, b>0 \theta>1 (-\Delta)_p^s denotes $-Laplacian operator 0<s<1 and 1<p<\infty N>sp \theta p<p_s^* p_s^* \frac{Np}{N-sp} critical Sobolev exponent, f superlinear continuous function subcritical growth V positive potential. Using penalization method Ljusternik-Schnirelmann theory, we study existence, multiplicity concentration of nontrivial solutions for small enough.